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ABSTRACT 

We use a stacking technique to measure the average Hi content of a volume-limited 
sample of 1871 AGN host galaxies from a parent sample of galaxies selected from 
the SDSS and GALEX imaging surveys with stellar masses greater than 1O 1O M0 and 
redshifts in the range 0.025 < z < 0.05. Hi data are available from the Arecibo Legacy 
Fast ALFA (ALFALFA) survey. In previous work, we found that the Hi gas fraction in 
galaxies correlates most strongly with the combination of optical/UV colour and stellar 
surface mass density. We therefore build a control sample of non-AGN matched to the 
AGN hosts in these two properties. We study trends in Hi gas mass fraction (Mhi/M*, 
where M* is the stellar mass) as a function of black hole accretion rate indicator 
L[OlIl]/Msjj- Wc find no significant difference in Hi content between AGN and control 
samples at all values of black hole accretion rate probed by the galaxies in our sample. 
This indicates that AGN do not influence the large-scale gaseous properties of galaxies 
in the local Universe. We have studied the variation in Hi mass fraction with black hole 
accretion rate in the blue and red galaxy populations. In the blue population, the Hi 
gas fraction is independent of accretion rate, indicating that accretion is not sensitive 
to the properties of the interstellar medium of the galaxy on large scales. However, 
in the red population accretion rate and gas fraction do correlate. The measured gas 
fractions in this population are not too different from the ones expected from a stellar 
mass loss origin, implying that the fuel supply in the red AGN population could be a 
mixture of mass loss from stars and gas present in disks. 

Key words: galaxies: evolution - radio lines: galaxies 
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INTRODUCTION 

Observations have shown that most low redshift mas- 
sive ga laxies con t ain a black hole (BH) in their 
bulges (jKormendvl 120041 ') and that in some cases, the 
black hole may be actively accreting. The mass of 
the black hole is strongly coupled with the mass 
and vel ocity dispersion of t h e bu lg e of its host 
galaxy (Kormcndy & Richstonc 1995; Magorrian ct al. 
19981: iGebhardt et all l200d : iFerrarese fcMerrittl |200d : 



Tremaine et all I2002T ). Whether or not energy output 



by accreting black holes can influence the evolution 
and properties of their host galaxies is still a mat- 
ter of debate. In some theoretical models of galaxy 



fabello@MPA-Garching.MPG.DE 



evolution, (radio) AGN are assumed to heat the ex- 
ternal gas reservoir in massive dark matter halos (eg . 
Bower et al1l2006l: ICattaneo et alj|2006l; ICroton et al.ll200d : 



De Lucia fc Blaizotl 120071 ; ISomerville et all |200S| ), thereby 
quenching the growth of the most massive galax- 
ies. It is believed that this process may solve the 
over-cooling problem in simulated galaxy clusters (eg . 
Churazov etafl l200ll: iRevnolds et all l200ll: lOuilis et al 



20021 : ISiiacki fc Sprin 



2001; Ruszkowski & Begelman 

2009 : iTevssier et all l201C| h Another favored evolutionary 
scenario is that the energy output by the AGN actu- 
ally drives gas out of the galaxy itself, depleting its in- 
terstellar medium, and shutting down sta r formation (e.g. 
iDi Matteo et al.ll2005l : iHopkins et al]|2006D . 

Observationally, the situat ion is still ve r y un- 
clear. In nearby ga l axies (e.g. iKauf fmann ct "all l2003al : 
iHeckman et all 20041 : ISchawinski et al.l 120071 ') and galax- 
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ies at intermediate redshifts (e.g. iNandra et al. 2007; 
ISilverman et al.ll2008l ). the host galaxies of AGN are found 
to occupy the "green" part of the colour - magnitude dia- 
gram, i.e. their stellar populations are intermediate in age 
between the blue, actively star-forming spiral galaxies, and 
the gas-poor, quiescent early-type population. One interpre- 
tation of these results is that there is an evolutionary se- 
quence driven by nuclear activity: i.e. AGN are triggered 
in galaxies on the blue sequence, gas in these systems is 
heated, and galaxies move onto the red sequ ence as star for- 
mation shuts down (Sc hawinski et all [20091. However, this 
interpretation is certainly not unique. At a fixed black hole 
mass, AGN tend to be found preferentially in galaxies with 
younger stellar populations (and hence more gas). The pres- 
ence of a gaseous fuel supply may simply be a necessary 
condition for both star formation and accretion onto the 
black hole (Kauffma nn et al"1l2003al ; iKauffmann fc Heckmanl 
2009). One way to ascertain whether AGN are responsible 
for terminating the growth of massive galaxies, is to estab- 
lish a direct link between the presence of an active nucleus 
and the gas content of such galaxies. 

There have only been a few systematic studies of the 
Hi conte nt of the host ga laxies of moderately large samples 
of AGN. iHo et al.l <|2008h conducted a survey of neutral gas 
in a representative sample of 154 moderate to high luminos- 
ity nearby Type I AGN. Their results appear to challenge 
the feedback scenario because AGN are found to be at least 
as ga s rich as quiescent objec ts of the same morphological 
type. [Hughes &: Cortesd |2009l ) carried out a study of Hi in a 
volume-limited sample of galaxies in different environments 
around the Virgo cluster, trying to disentangle between the 
effect of AGN and environment on star formation. In their 
sample, the incidence of AGN peaks in the "green valley", 
but galaxies with quenched star formation are found to be 
mainly HI deficient galaxies located in high-density envi- 
ronments. No direct connection between AGN activity and 
quenching of star formation was observed, at odds with the 
previous cited studies. 

In this paper we measure the average Hi content of 
a sample of ~2000 AG N host galax i es, us ing the stacking 
techni que developed in |Fabello et ail l|2010l . hereafter Paper 
I). In ICatinella et al.l and Paper I we showed that 

the gas content of galaxies with stellar masses greater than 
1O 1O M0 is, to first order, most tightly correlated with op- 
tical/UV (NUV— r) colour and stellar surface mass density. 
We build control samples that are very closely matched to 
our AGN sample in the (NUV-r)//i, plane, and we look for 
differences in Hi mass fraction Mhi/M* between AGN and 
control galaxies. In section 1 we describe the sample selec- 
tion and the data we use; in section 2 we compare in bins 
of nuclear properties the gas content of AGN and inactive 
counterparts. 



1 SAMPLES 

Our targets are drawn from the volume-limited, homoge- 
neous sample A defined in Paper I. To summarize, this sam- 
ple consists of 4726 galaxies with stellar masses greater than 
1O 1O M and redshifts in the range 0.025 < z < 0.05 selected 
from the SDSS main spectroscopic sam ple and the projecte d 
GALEX Medium Imaging Survey (MIS; lMartin et al.ll2005h . 



HI data are availab le from the Arecibo L egacy Fast ALFA 
burvey (ALFALFA; iGiovanelli et~ali r2005). Among our tar- 
gets, 23% have been detected by ALFALFA. We refer the 
reader to Paper I for further details on the sample selec- 
tion and for a description of the stacking method used to 
evaluate average Hi gas fractions. 

From this main sample, we extract a sample of AGN 
(5 11. 2p , and a control sample of galaxies where there is no 
evidence of accretion onto the central black hole f§ ll.3p . In 
the next section, we introduce the parameters we use for 
both the analysis and the sample selection. 

1.1 Galaxy Parameters 

The optical/UV parameters we use are derived from SDSS 
spectru m measurement Jj] a nd reprocessed GALEX UV pho- 
tometry rtWang et al.ll2010l . Paper I). From there we use: the 
NUV— r colour measured for the whole galaxy (corrected 
for Galactic extinction only); the stellar mass surface den- 
sity /x*=M*/(27r_R§ z ), where Rso,z is the Petrosian radius 
containing 50% of the flux in z-band in units of kpc; the 
central velocity dispersion a. In addition we use the 4000A 
break measured from the SDSS 3 arcsec-diameter fiber spec - 
trum, D n (4000). As discussed in IKauffmann et al.l (|2003bh . 
-D„(4000) can serve as a measure of the age of the stellar 
population of galaxies, although there is a secondary depen- 
dence on metallicity. The galaxies in our sample are located 
at low redshifts (0.025 < z < 0.05), and in 72% of our sam- 
ple the 3 arcsec fiber covers less than 20% of the Petrosian 
radius enclosing 90% of the total light of the galaxy. The 
4000A break values should thus be regarded as a measure 
of stellar population age for the inner nuclear regions of the 
galaxies, while the NUV— r colour is a proxy for the specific 
star formation rate SFR (sSFR, i.e. SFR per unit stellar 
mass) . 

Finally, we extract the following emission fluxes: Ha, 
H/3, [Nil], [Om] (in units of 10" 17 erg-s _1 cm~ 2 ), which we 
use to classify the nuclear activity of our gal axies by means 
of the diagnostic diagram first defined by iBaldwin et al.l 
l|l98ll. hereafter BPT). Because the fiber covers the inner 
region of the galaxy, the AGN properties should not be sig- 
nificantly contamin ated by outer star forming regions (see 
iKewlev et aTll2006l ). 

As discussed in IKauffmann et alj (|2003aT ) ■ the luminos- 
ity of the [Om]A5007 line should be a reasonably reliable 
tracer of black hole accretion rate. Because it is a high ex- 
citation line, it is less contaminated by emission from Hn 
regions than other lines such as Ha. We correct our line 
fluxes for dust ex tinction usi n g the B aimer decrement, fol- 
lowing eq. (3) in IWild et alj 1 2007h . and then convert the 
[Om] fluxes into luminosities. 

Finally, we evaluate an Eddington ratio for the AGN, 
defined as: 

L bol L bo i L[OIII] 

(X — oc — 

Ledd Mbh & 

where Lbol is the bolometric luminosity, and Mbh the cen- 
tral black hole mass. The Eddington luminosity Ledd scales 

1 See http://www.mpa-garching.mpg.de/SDSS/DR7/ and 
http: / / cas.sdss.org. / dr7/en/tools/search/sql.asp 
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Figure 1. BPT diagram for galaxies in our sample with S/N>3 
for the four emission lines Ha, H/3, [Oin] and [Nil]. The solid curve 
we use to demarcate the boundar y between AGN and "no rmal" 
star- forming (SF) galaxies is from lKauffmann et al. | j2003aT ); the 
dash ed line shows the d emarcation boundary for "pure" AGN 
from lKewlev et al.l l l200ll) . Objects are colour-coded according to 
their nuclear properties as labeled in the diagram. 



with Mbh, which we estimate using the stellar velocity dis- 
persion a measured from the fiber spectrum: M bh oc a 4 . 
Lboi scales with L[Oiii] fsee lHeckman et al.| [2004). so to first 
order the ratio L[Oiii]/<t 4 can be used as a proxy for the 
accretion rate. For a precise est imate of the Eddingto n pa- 
rameter we refer the reader t o iHeckman et all l|2004 ) and 
iKauffmann fc Heckmanl (|2009h . 

1.2 The AGN sample 

We extract a sample of AGN from sample A using the 
BPT diagnostic diagram, which allows to distinguish Type 
2 AGN from star-forming galaxies. We consider the ratio 
[Oili]/H/3 versus [Nil] /Ha for all galaxies with signal-to- 
noise ratio S/N > 3 in al l four line measurements. Following 
IKauffmann et al.1 i|2003aT ). AGN are defined as those galaxies 
with 

log ([OIH]/H#) > 0.61/{log([NII]/Ha) - 0.05} + 1.3, 

We obtain a sample of 1871 active galaxies, which can be 
further subdivided into Seyferts, LINERs (which make up 
the bulk of our sample) or "transition" as indicated in Figure 

m 



We note that iKewlev et all (|200ll ) suggested a more 
stringent cut to separate galaxies where the nuclear emis- 
sion is almost completely excited by emission from the gas 
accreting onto the black holes, rather than from star-forming 
regions (shown as a dashed line in Figure [lj. With this cut, 
we obtain a subsample of 912 "pure" AGN, which is too 
small for our stacking analysis. T his is why we only consider 
the demarcation line suggested bv lKauffmann et al.l (|2003al ) 
in this paper. 

In the top panel of Figure [2j we plot the AGN in the 
colour - stellar ma ss surface dens i ty pla ne. The "pure" AGN 
defined using the IKewlev et al.l (|200ll ) cut (dashed line in 
Figure [TJ are over-plotted as magenta points. As can be 
seen, both samples span the same range of parameter space. 
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Figure 2. Galaxies in our samples are plotted in the 
(NUV-r)/)i* plane. Top: Black points show the AGN sample. We 
overplot as magenta points galaxies that lie above the stronger cut 
defined by IKewlev et all ll200ll) . Bottom: All control pool galax- 
ies, made by "inactive" galaxies with S/N([Om])< 3 plus star- 
forming objects. Gray dots are the galaxies that are discarded 
when matching to the AGN sample (i.e. mainly very blue, low 
stellar surface density objects). Black and coloured dots show the 
final control sample extracted. In particular, galaxies for which all 
four emission lines have S/N < 3 are indicated as red points, while 
the blue points show galaxies that are classified as star-forming. 



1.3 The control sample 

Our control sample is drawn from a combined sample of "in- 
active" galaxies and star-forming systems. Inactive galaxies 
are those where the S/N of the [Oin] line is less than 3. 
We exclude 474 galaxies where the S/N of the [Oin] line is 
greater but other lines lie below this S/N threshold, because 
these galaxies cannot be accurately classified. By applying 
a cut in the [Om] line only, our control sample may be con- 
taminated by weak LINERs. The contamination is anyway 
negligible. For our targets, less than 5% of the galaxies with 
S/N([Om]) < 3 have log([NII]/Ha) > log(0.6) (vertical dash 
dotted line in Figure [T}. 

Inactive galaxies, together with galaxies classified as star- 
forming on the BPT diagram, constitute our pool of 2377 
galaxies for which the nuclear emission is not dominated by 
an AGN. The colour-^* plane for those galaxies is shown in 
the bottom panel of Figure [2] all points. Two distinct popu- 
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lations are clear in this diagram: a red sequence of galaxies 
with weak or absent emission lines, and a blue cloud con- 
sisting of star-forming galaxies which extends until lower 
with respect to the active sample. 

We stack AGN in bins of the accretion rate proxy 
L[Oin]/<7 4 . For each bin, we want to compare the aver- 
age Hi mass fraction with that obtained for control galaxies 
matched in stellar mass surface density and NUV— r colour. 
For each AGN we then select a galaxy from the control pool 
by searching for its closest neighbor in the plane defined by 
NUV— r and /i*. We match objects in order of decreasing 

and we do not allow a galaxy to enter the control sam- 
ple more than once. We do not have enough galaxies to build 
larger control samples, but with this method we are able to 
reproduce the original and NUV— r distributions of the 
AGN sample. The gray dots in Figure [2] (bottom panel) are 
the objects discarded because they do not match any AGN. 
The black and coloured points show the final control sam- 
ple. In particular, blue dots are star forming galaxies, while 
galaxies for which none of the four emission lines Ha, H/3, 
[Oin] and [Nil] are detected with S/N > 3 are overplotted 
as red dots. 



2 THE HI GAS FRACTIONS OF AGN HOSTS 

We have stacked our AGN sample in bins of [Oin] line lu- 
minosity and in bins of the black hole accretion rate proxy 
L[Oiii]/<t 4 . In Figure [3] we plot the average gas fraction 
Mhi/M* as a function of both parameters for both the AGN 
sample (red) and for the control galaxies (black). Errors are 
evaluated using the bootstrap method for a total of 200 rep- 
etitions: every time we discard 20% of the sample and eval- 
uate a new gas fraction; the error is the dispersion on the 
repeated measures. 

As can be seen, the average Hi mass fraction increases 
weakly as a function of [Om], from about 9% for the weakest 
AGN to ~14% for the most luminous systems. This nearly 
flat trend reflects the fact that AGN of different [Om] lu- 
minosities are located in roughly the same region of the /x*- 
(NUV— r) plane. The dependence of Hi mass fraction on 
black hole accretion rate is considerably stronger. Objects 
with the lowest measurable accretion rates are very gas poor 
(Mhi/M*~3%), and the average gas fraction increases as a 
function of L[Oiii]/<7 4 up to a maximum value of ~14%. We 
likewise find that galaxies with higher black hole accretion 
rates tend to be bluer and have lower values of /i*. We note 
that the trend in Hi mass fraction as a function of black 
hole accretion rate is considerably weaker than as a func- 
tion of the optical/UV colour of the host galaxy. In Paper 
I, we found that the average Hi mass fractions change from 
~55% for the bluest objects to ~2% for the reddest ones. 

The next step is to assess whether AGN differ from con- 
trol galaxies in terms of mean Hi gas mass fraction. Results 
for the control sample are plotted as black stars in Figure 
[3] at the x-axis value of the corresponding AGN bin. Errors 
are evaluated as explained previously. As can be seen, there 
is no difference between the neutral gas content for the two 
sam ples. 

Kauffmann fc Heckmanl (|2009h showed that there ap- 
pear to be two distinct regimes of black hole growth in galax- 
ies in the local Universe. The first is associated with galaxies 



with significant ongoing star formation; in this regime the 
distribution of accretion rates shows little dependence on the 
central stellar population of the galaxy. The second regime 
is associated with "passive" galaxies, and is characterized by 
a power-law distribution of accretion rates. There, the ac- 
cretion rate does depend strongly on the age of the central 
stellar population in the galaxy. 

We now take a look at Hi gas mass fraction trends for 
blue and red galaxy sub-populations. We define red sequence 
galaxies to have NUV— r > 4.5, and blue sequence galaxies 
to have NUV— r < 4.5. This cut is based on the split in the 
bimodal colour distribution of the control sample, visible 
also in Figure [2] (lower panel). Results are shown in Figure [4] 
(left panel): blue stars are blue cloud AGN, black stars are 
the corresponding control galaxies; red triangles represent 
the red sequence AGN, while black triangles are used for 
control galaxies. Symbols are empty when the measure is a 
non-detection (upper limit). The average Hi gas fraction re- 
mains constant as a function of the accretion rate parameter 
L[Oiii]/<t 4 for the blue population. The increase in Hi gas 
fraction as a function of L[Om]/cr 4 seen in Figure is in fact 
driven by galaxies on the red sequence. We note that we had 
to increase the bin size for these objects in order to recover 
sufficient signal in our stacks. Our control samples of red and 
blue galaxies exhibit the same trends with L[Oiii]/<7 4 as the 
AGN, but the gas mass fractions are slightly higher for the 
blue control objects, and lower for the red control objects. 
(In fact, we do not detect Hi in any of the red control galaxy 
stacks). 

The NUV— r colours we use are not corrected for dust 
extinction. It is conceivable that AGN may be found in 
galaxies with more dust than average, so we should ascertain 
whether using uncorrected colours to create matched control 
samples will bias our results. Reliable estimates of the spe- 
cific star formation rates of our galaxies may be obtained 
by fitting spectral energy distribution models to the 7-band 
UV and SDSS data, and calibrating the attenuation as func- 
tion of NUV— r colour using a reference sample with direct 
measurements of the UV through far-IR spectral energy dis- 
tribution from a combination of GALEX, SDSS and Spitzer 
data. Our methodology is explained in detail in Saintonge 
et al. (2011, MNRAS, submitted). 

We repeat our stacking exercise, this time matching 
AGN hosts to control galaxies using /i* and specific star 
formation rate rather than [i+ and NUV-r colour. The 
new results are shown in the right panel of Figure [4] Red 
symbols denote objects with low specific star formation 
rates (logSFR/M* < -11.0), while blue symbols are for more 
strongly star-forming galaxies (log SFR/M* > -11.0). Once 
again, this division has been made based on the bimodal dis- 
tribution of sSFR in our sample. We see the same trends as 
before, with no significant difference in HI content between 
the blue AGN hosts and control galaxies. The same conclu- 
sion holds for red AGN hosts, although we note that the two 
highest accretion rate bins yield non-detections. This is be- 
cause many of these galaxies shift over to the blue sequence 
after dust corrections are made. 

We note that the amount of accretion onto the cen- 
tral supermassive black hole in a galaxy is more likely to 
be sensitive to the properties of the interstellar medium 
in the central parts of the galaxy rather then globally. We 
now investigate what happens if we split our AGN sample 
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Figure 3. Hi gas fraction as a function of L[Om] (left panel) and L[Oiii]/<t 4 (right panel). Filled red stars are the values obtained for 
the AGN hosts, while the black stars are the values obtained for the corresponding control galaxies. The number of galaxies included in 
each stack is reported at the bottom of the plot at the position of each bin. 
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Figure 4. Hi gas fraction as a function of L[Oiii]/<7 4 for galaxies and AGN hosts, split into two bins of star formation properties. Coloured 
symbols show results for galaxies with active galactic nuclei, while black symbols show results obtained for the corresponding control 
samples. The empty symbols denote measurements that are upper limits. The number of objects stacked in each bin is reported. Left 
panel: the trend with gas fraction is shown separately for blue galaxies (defined as NUV— r<4.5), and red sequence ones (NUV— r>4.5). 
Right panel: galaxies are matched and divided into two bins of specific star formation rate; blue symbols are for star forming targets, 
red for quiescent ones. 



into "red" and "blue" subsamples using D n (4000), which 
is measured within the 3 arcsec diameter SDSS fiber aper- 
ture, and remeasure gas fractions for these two new samples, 
plus corresponding control galaxies. In this case, blue galax- 
ies have D„(4000)<1.6 and red galaxies have D n (4000)>1.7 
(this roughly separates the peaks of the bimodal distribu- 
tion of 4000A break strengths) . In Figure [5] we show our 
results. AGN hosts exhibit no significant differences in their 
gas content with respect to non-AGN in the same D n (4000) 
range matched by global colour and fi±. Also the differences 
in the red population seen in Figure 2] has now vanished, 
given that now we detect the same signal for both AGN 
and control galaxies. We note that the two bins of D„ (4000) 
clearly separate in terms of the trend in gas fraction as a 
function of the Edding ton parameter L[Olll1/q- 4 , cons istent 
with the suggestion of iKauffmann fc Heckmar] |2009l ) that 



accretion onto the black hole is supply-limited in the red 
galaxy population, but not in the blue one. 



3 SUMMARY & DISCUSSION 

We have used the stacking technique developed in a compan- 
ion paper (Paper I) to look at possible interplay between 
nuclear activity and gas content in massive nearby galax- 
ies. Starting from a volume-limited, homogeneous sample of 
massive galaxies we have selected a sample of AGN, and 
a control sample of inactive galaxies matched by and 
NUV— r. We then stacked their Hi data in bins of [Om] line 
luminosity scaled by black hole mass to assess whether the 
gas content of AGN hosts differs from that of the control 
galaxies. 

Our main result is that we do not find any significant 
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Figure 5. Hi gas fraction as a function of L[Oiii]/<t 4 , for galaxies 
divided into two subsamples according to 4000A break strength: 
red triangles are galaxies with D n (4000) > 1.7 and blue stars with 
D n (4000) < 1.6. Symbols are as in Figure [5] 



difference between the Hi gas content of AGN hosts and that 
of control galaxies. This conclusion holds at all values of the 
Eddington parameter proxy L[Oiii] /a 4 probed by the galax- 
ies in our sample. 

If converted into a proper Eddington ratio, our sample would 
cover the range log(L/LEdd) ~-4 to -2. We are studying the 
relatively low accretion rate regime, and that could be a rea- 
son why we do not see any sign of AGN feedback. We would 
need data over a larger volume to stack a significant num- 
ber of AGN with very high L/IjEdd- Still, up to the highest 
accretion rates we have studied, the gas fraction shows no 
hint of correlation with either L[Oiii]/<t 4 (blue population) 
or AGN presence. 

These results are also in agreement with iHo" et al.l 

(2008). As discussed in their paper, the lack of any differ- 
ence may be due to the fact that the Hi traces gas on very 
large scales, but accreting black holes only heat gas close to 
the centers of the galaxies. 

In some simulations of galaxy formation, AGN feed- 
back is an extremely violent process that can drive out much 
of the interstellar medium from galaxies, truncating ongo- 
ing star formation and reddening their stellar populations 
i|Di Matteo et al.ll2005l ). In this scenario, AGN hosts would 
be expected to be gas deficient. Catastrophic global feed- 
back processes are clearly not occurring in the majority of 
AGN in the local Universe. We cannot yet ascertain that 
AGN feedback is important in the centers of galaxies. In 
future work, we plan to examine trends in molecular gas 
fractions in AGN hosts using data from the COLD GASS 
survey (Saintonge et al 2011). The molecular gas is generally 
more centrally concentrated that the atomic gas and may be 
more sensitive to feedback processes occurring in the bulge 
of the galaxy. 

Another interesting result of our analysis is the scaling 
of the mean Hi gas fraction of AGN host galaxies with ac- 
cretion rate onto the black hole, as traced by L[Olll] /a 4 . In 
the blue galaxy population, the accretion rate onto the black 
hole and the Hi gas fraction of the galaxy are independent. 
However, in the red galaxy population, these two q uantities 
do track each other. iKauffmann fc Heckmanl (|2009l ) hypoth- 



esized that in blue galaxies, the black hole regulates its own 
growth at a rate that does not further depend on the proper- 
ties of the interstellar medium. In red galaxies, they deduced 
that the decrease in the accretion rate onto black holes in 
old galaxies was consistent with population synthesis model 
predictions of the decline in stellar mass loss rates as a func- 
tion of mean stellar age. 

Because we now have measured Hi mass fractions for 
the red population, we can check this in more detail. The 
average global (NUV— r) /{g — r) colours of the host galaxies 
vary from 4.56/0.74 for the stack with the highest accretion 
rate in Figure[S]to 5.54/0.80 for the stack with the lowest ac - 
cretion rate. As described in IKauffmann fc Heckmanl (120091) 
we use population syn thesis models (e.g-. lBruzual fc Charlotl 
l2003l;lMarastonll2005l ) to calculate mass loss rates as a func- 
tion of colour, assuming simple exponential declining star 
formation histories for our galaxies. If the Hi content of red 
sequence AGN came only from mass loss, we would expect a 
0.3 dex change in the average log (Mhi/M*) between AGN 
stacks with the lowest and highest accretion rates in our 
sample. Figure [S] shows a 0.5 dex change in this quantity, 
which is somewhat larger than predicted. However, we note 
that the average stellar surface mass densities of the host 
galaxies in the stacks decrease by a factor of 3 between AGN 
with the highest and lowest accretion rates (likewise the con- 
centration index decreases from 3.15 to 2.79), so structural 
properties are not constant across the population. This may 
imply that the fuel supply in the red AGN population is in 
fact a mixture of mass loss from stars and gas present in 
disks. 
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